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ABSTRACT 

In this paper we combine high resolution N-body simulations with a semi analytical model of galaxy 
formation to study the effects of a possible Warm Dark Matter (WDM) component on the observable 
properties of galaxies. We compare three WDM models with a dark matter mass of 0.5, 0.75 and 
2.0 keV, with the standard Cold Dark Matter case. For a fixed set of parameters describing the 
baryonic physics the WDM models predict less galaxies at low (stellar) masses, as expected due to the 
suppression of power on small scales, while no substantial difference is found at the high mass end. 
However these differences in the stellar mass function, vanish when different set of parameters are 
used to describe the (largely unknown) galaxy formation processes. We show that is possible to break 
this degeneracy between DM properties and the parameterization of baryonic physics by combining 
observations on the stellar mass function with the Tully-Fisher relation (the relation between stellar 
mass and the rotation velocity at large galactic radii as probed by resolved HI rotation curves). WDM 
models with a too warm candidate (m„ < 0.75 keV) cannot simultaneously reproduce the stellar mass 
function and the Tully-Fisher relation. We conclude that accurate measurements of the galaxy stellar 
mass function and the link between galaxies and dark matter haloes down to the very low-mass end 
can give very tight constraints on the nature of DM candidates. 

Subject headings: methods: analytical - galaxies: mass function - galaxies: formation - cosmology: 
theory - dark matter - large-scales structure of Universe 



1. INTRODUCTION 

The Cold Dark Matter (CDM) paradigm successfully 
describes the formation of large-scale structure in the 
Universe (e.g., Springel et al. 2006). On small scales, 
the CDM model faces inconsistencies with observations. 
For example, CDM model over-predicts the number of 
satellites in the Milky Way (e.g. Klypin et al. 1999), 
predicts cuspy density profile for the dark matter halo 
(e.g. Moore et al. 1999), and a too large kinematics for 
massive satellites (Boylan-Kolchin et al. 2011). 

These small-scale problems in CDM all arise from the 
comparison of pure N-body simulations with observa- 
tional data. Recently it has been shown that some of 
those problems could be alleviate by the inclusion of 
baryonic physics, due to its back-reaction on the proper- 
ties of dark matter haloes. For example, UV background, 
rcionization and supernova feedback can act together to 
suppress the formation of dwarf satellites in the Milky 
Way (e.g. Bullock et al 2001, Maccio et al. 2010), the 
inner halo density profile becomes more flat when star 
formation and feedback are included (e.g. Mashchenko et 
al. 2008, Pontzen& Governato 2012, Maccio etal. 2012). 
However, these baryonic effects are difficult to model and 
still strongly rely on simple parameterization for complex 
physical processes (the so called sub-grid physics) 

Another possible solution is to suppress the excess of 
power at small scales in the Dark Matter distribution by 
changing the properties of the dark matter candidate. 
Among many of these different approaches, the Warm 
Dark Matter (WDM) model is the most intriguing as it 
naturally preserves the large-scale structure of the CDM 
model (Alam et al.2002; Zentner & Bullock 2002). Many 
authors have shown that the WDM model can relax some 



tensions between observations and theoretical predic- 
tions using pure N-body simulations (Colin et al. 2000, 
Jing 2001; Bode et al. 2001, Gotz & Sommer-Larsen 
2002, Knebe et al. 2008, Colin et al. 2008, Tikhonov et 
al 2009, Lovell et al.2012, Maccio et al. 2012b). 

One possible WDM candidate is a sterile neutrino 
which exhibit a significant primordial velocity distribu- 
tion and thus damp primordial inhomogeneities on small 
scales (e.g. Hansen et al. 2002, Abazajian & Koushiap- 
pas 2006, Boyarsky et al. 2009). Limits on the mass of 
dark matter particles can be obtained from several as- 
trophysical observations: one of the most powerful tool 
for constraining the matter power spectrum are Lyman-a 
forest observations (neutral hydrogen absorption in the 
spectra of distant quasars, Narayanan et al. 2000, Vicl 
et al. 2005, Seljack et al. 2006). Lyman-a observations 
allows the possibility to studying the power spectrum 
down to small scales and over a large range of redshifts. 
Current observations set a lower limit of m„ « 1 keV 
(Viel et al. 2008). 

These limits have been confirmed by other observations 
based on different methods like QSO lensing (Miranda & 
Maccio 2007), luminosity function of high redshift QSOs 
(Song & Lee 2009), dwarf galaxies in the Local Volume 
(Zavala et al. 2009), the size of (min)-voids around the 
Local Group (Tikhovov et al. 2009) and Milky Way satel- 
lites (Maccio & Fontanot 2010). 

Up to now, most studies based on numerical simula- 
tions have focused their attention on the internal struc- 
ture and kinematics of nearby galaxies in the Warm Dark 
Matter model, while very few have explored the effect 
that a warm candidate will have on the properties of the 
general galaxy population, such as the luminosity func- 



2 



tion, the stellar to halo mass function and so on. 

A notable exception is the recent paper by Menci 
et al. (2012), although not directly based on N-body 
simulations in WDM. In their study the authors com- 
bined a Semi Analytical Model for galaxy formation with 
the Extended Press and Shechter formalism (e.g. Cole 
et al. 2000, Lacey & Cole 1993) to obtain the properties 
of galaxy in a WDM model. They found that the WDM 
model can produce a more flat slope for the faint end 
of the luminosity function in better agreement with the 
observations. They also argue that other galaxy proper- 
ties better match observations when a warm candidate 
is used instead than a cold one. While these effects are 
certainly true, the adoption of a Warm Dark Matter like 
spectrum has also strong influence on the ratio between 
dark matter and stellar mass in galaxies. This ratio that 
can be constrained by several dynamical mass estimators 
such as the Tully-Fisher relation (Tully & Fisher 1977), 
satellites kinematics, and weak gravitational lensing. 

In this paper we combine high resolution N-body sim- 
ulations for the standard CDM scenario and three WDM 
candidates = 2,0.75,0.5 keV, with the Semi Analyt- 
ical Models first presented in Kang et al. (2005), and 
then successively expanded and improved in Kang & van 
den Bosch (2008) and Kang et al. (2012). Our aim is 
to study the impact of different WDM models on galaxy 
properties and use different observations to try to disen- 
tangle the effects of a warm candidate from the different 
parameterization of the baryonic physics adopted in the 
SAM. The final goal is to constrain the properties of the 
dark matter candidate. 

The paper is organized as follows: We introduce the 
simulations and model in §2, show the model predictions 
in comparison to observations in §3, and conclude in §4. 

2. N-BODY SIMULATIONS AND GALAXY FORMATION 
MODEL 

In the WDM model, the velocity dispersion of dark 
matter particles defines a characteristic free-streaming 
scale, below which the perturbation will be suppressed. 
The free-streaming scale can be described as (Bode 
et al. 2001), 

fn m h 2 \ 1/3 fm v \~ 4 /3 
^ = °- n (w) (keV) MPC W 

The lighter the particle mass is, the larger the scale 
below which perturbations will be suppressed. It also 
possible to define a free-streaming mass scale as, Mf s = 
4/37r ( o c (rf s /2) 3 (where p c is the critical density of the 
Universe), and the initial density perturbations are prac- 
tically erased on masses below Mf s . The impact of the 
free-streaming on the power spectrum can be described 
using a modified transfer function as suggested in Bond 
et al. (2001): 

T (k) = (TOTrJ = C 1 + («fc) 2A T 5Al - (2) 

Where -fLi^ DM and are the linear power spectrum 

in the WDM and CDM models. Viel et al.(2005) using 
a Boltzmann simulation code found that /i = 1.2, and 
a can be directly related to the mass of the warm dark 
matter particles and the matter density parameter in the 



following way: 

a = 0.049 f^r 111 ( ikV U ( ±) fc -i M pc. 
VlkeV/ V.0-25 / \0.7J v 

We used the above equations, to derive the initial 
power spectrum for our WDM models, and take it as 
input for generating the initial conditions of our N-body 
simulation. 

The cosmological parameters are set in agreement with 
those derived from the WMAP7 data release (Komatsu 
et al. 2011), namely: ft A =0.73, Q m =0.27, £4=0.044, 
h = 0.7 and as = 0.81. The simulations have been run 
using the GADGET-2 code (Springel 2005) in a box of 
L = 200 Mpc/h using 1024 3 particles. This sets our 
mass resolution to 5.5 x 10 s /i _1 M Q and the softening to 
4.5 /i kpc. We run three different WDM models for a 
thermal candidate with particle mass of 0.5, 0.75 and 2.0 
keV respectively. Finally, readers are refereed to Kang 
et al. (2005) for details on the halo finder and how we 
construct merger trees from the simulations. 

The Semi Analytical Model (SAM) of galaxy formation 
we employ in this study is based on the model of Kang 
et al. (2005, 2008). This SAM self-consistently models 
the physical processes governing stellar mass evolution, 
such as gas cooling, star formation, supernova and active 
galactic nucleus (AGN) feedback. The latest incarna- 
tion of this model had been recently presented in Kang 
et al. (2012) and we refer the reader to this paper for 
more details. As in Kang et al. (2012), we slightly mod- 
ified the cooling rate of low-mass haloes as, 

Aicool = fc X mhotAdyn- (4) 

where f c is the cooling factor. As described in Kang 
et al. (2012), this formula can be interpreted as an effec- 
tive rate for gas cooling, that takes into account the gas 
outflow due to reheating by supernova feedback. 

Kang et al. (2012) has shown that a low value for f c is 
needed to better reproduce the slope of the low-mass end 
of the galaxy stellar mass function in the CDM model. 
We will show later in §3 that for the WDM models, due 
to lower power on small scales, this cooling factor f c has 
to be increased to match the local stellar mass function. 

3. RESULTS 

In the WDM model, the formation of haloes with mass 
below the free-streaming scales is suppressed. For the 
model with m„ = 0.5 keV, the free-streaming mass is 
about 10 9 M Q . Actually also the formation of halo more 
massive than this cut-off can also be affected up to a 
mass of the order of 2000 x Mf s (e.g. Schneider et al. 
2012). 

In the left panel of Fig. [T]we show the dark matter halo 
mass functions for our different models. Note that this 
mass function is slightly different from the usual mass 
function that normally only includes Friends Of Friends 
(FOF) haloes (i.e. central haloes). In our mass func- 
tion we also include satellite halos: their mass is defined 
as the mass at the accretion time M acc , in other words 
the mass the satellite had before becoming a satellite of 
a more massive halo. It has been shown that this ac- 
cretion mass is more closely related to the final stellar 
mass of the satellite (e.g., Vale & Ostriker 2004, Moster 
et al. 2010). For central halos, the mass is simply defined 
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Fig. 1. — The halo and stellar mass functions. Left panel: the halo mass (M aC c) functions from the CDM and WDM simulations . Middle 
and right panels: galaxy stellar mass function at z = and z = 1.5, lines are from our model with / c = 0.15 (see text). The z = data 
points are from Li & White (2009) and Bell et al. (2003), and high-z data points are from Marchesini et al. (2009), Fontana ct al. (2006). 



as the virial mass at the current time. Such a mass func- 
tion which also includes satellite galaxies can be easily 
linked to the predicted galaxy mass functions, which is 
one of the observations we want to study in Warm Dark 
Matter. 

For our more extreme WDM model with m„ = 0.5 
kcV, the halo number density is lower by roughly 60% at 
Mdm = 10 1 M©. This confirms that the free-streaming 
in WDM has non-negligible effects on the formation of 
haloes with mass much larger than Mf s . On the other 
hand the model with m v = 2 keV, has a cutoff mass of 



10 M Pl which is well below the mass resolution 



of our simulations, as a consequence we do no see any 
differences with respect to the CDM model on the whole 
mass range probed by our runs. The warmest model 
also presents an upturn in the halo mass function around 
Mdm < 10 11,5 M Q . This upturn is due to the artificial 
clumping effect in WDM simulations from the limited 
resolution (Wang & White 2007; Schneider et al. 2012), 
which indicates that low mass haloes are possibly over- 
predicted in our simulations. We will later briefly discuss 
the implication of this effect on the prediction of galaxy 
properties. Finally all model agree reasonably well for 
masses higher than ~ 10 12 5 M Q , so we do not expect 
large differences for the properties of galaxies living in 
massive haloes. 

By applying the Semi Analytical Model to the N-body 
simulations it is possible to compute properties of the 
galaxies harbored in the dark matter haloes in the various 
WDM and CDM models. The middle and right panels 
in Fig. Q] show the stellar mass functions at z = and 
z = 1.5 respectively. In this first attempt we apply the 
same SAM (with the same parameters) to all the models. 

At z = the WDM model with lower m v produces 
too few galaxies around M sta r = 10 10 M Q with respect to 
current observational data; this is a direct consequence of 
the lower normalization to the halos mass function below 
Mdm = 10 n - 5 M Q in WDM with respect to CDM. The 
model for 0.75 keV gives a stellar mass function between 
the previous model and the CDM prediction, while the 
model with m„ = 2 keV produces almost identical results 
to the CDM model (as expected from the results on the 
halo mass function) and we will not show its results in 



the rest of the paper. The fact that WDM models seem 
to produce too few galaxies at low stellar masses with 
respect to observations docs not come as a surprise. The 
parameters in the SAM we use have been tuned to repro- 
duce the stellar mass function at z = in the (standard) 
CDM scenario, so any models that predict less low mass 
dark matter haloes will, by construction, underpredict 
the number of galaxies at low stellar mass. 

Results are somewhat inverted at higher rcdshift (z = 
1.5), as shown in the right panel of Fig [1] The CDM 
model ovcrprcdicts the number of galaxies with a stellar 
mass below 10 10 ' 5 M Q . This is a well know problem of 
the current modeling of galaxy formation (e.g., Fontanot 
ct al. 2010, sec Weinmann ct al. 2012 for a thorough de- 
scription of this issue). In this case having a warm dark 
matter component naturally helps in reducing the num- 
ber of low mass haloes, and consequently the number 
of low (stellar) mass galaxies, bringing SAM results in 
better agreement with observations. But of course this 
"success" comes at the sacrifice of local (z = 0) agree- 
ment. 

Recently Menci et al. (2012) used this suppression of 
stellar mass function at high rcdshift in warm dark mat- 
ter model as a hint that WDM may help in reducing the 
tension between data and theoretical prediction at high 
rcdshift. In their case the adopted SAM was also over- 
predicting the luminosity function at z = for the CDM 
model, where on the contrary our SAM+CDM model 
provides a good fit to the local data. So having a WDM 
component in the case of Menci et al. (2012) was helping 
on both high and low redshift. 

What remains to be clarified is whether the success of 
WDM claimed by these authors is truly due to the WDM 
component or simply reflects a non-optimal parameter- 
ization of the physics of galaxy formation implemented 
in their SAM (that could partially be the explanation of 
the failure of the CDM scenario in their model at both 
high and low luminosities, at z = and z = 1.5). 

In the same paper Menci and collaborators also found 
an effect of WDM on the high luminosity tail of the lumi- 
nosity function. Our results do now show such an effect 
- all our models give the exact same prediction at high 
stellar masses. This result is somehow expected give the 
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Fig. 2. — The galaxy stellar mass functions at z = (left) and z = 1.5 (right). As in Fig. [T] the points show the observations and the 
lines show the models. For the WDM models the cooling factor is treated as a free parameter. WDM models with m„ = 0.75, 0.5 keV, 
have f c = 0.35, 1.0 respectively. 



similarity of the halo mass function at high masses (see 
left panel in Fig. [lj. Mcnci et al. claimed that this dif- 
ference was due to the different accretion rates of satel- 
lites in massive haloes. SAMs have been show to have 
problems in modelling the properties of satellites in high 
mass haloes (e.g., Weinmann et al. 2010, Liu et al. 2010 
and references therein). This again raises the question 
whether the effect seen by Menci and collaborators (but 
not in this work) is a true effect due to warm dark mat- 
ter, or simply reflects the specific way in which baryonic 
physics is implemented in their semi analytical model. 
We will address this issue of degeneracy between SAM 
parameters and DM properties in the next section. 

3.1. Degeneracy between SAM parameters and WDM 

The failure of WDM models in matching the redshift 
zero stellar mass function does not imply that WDM 
models are somehow "wrong" ; the reason for this resides 
in the large uncertainty of baryonic physics parameteri- 
zation in the SAMs. Obviously the same argument ap- 
plies to the "success" of WDM models in reproducing 
the observational data at z = 1.5. 

The middle panel of Fig. Q] suggests that the stellar 
mass in low-mass halo is too low in WDM models. One 
possible solution would be to increase the star formation 
efficiency or decrease the effect the supernova feedback 
in the model, but this affect the high mass tail of the 
stellar mass function, where all models perform equally 
well. 

Instead we decided to turn to the effective cooling fac- 
tor given in eq. 21 since this factor is introduced only for 
low-mass haloes. To get a better match to the local stel- 
lar mass function, we set f c = 1.0 and 0.3 in the WDM 
models with m„ = 0.5 and 0.75 keV, respectively. 

Fig-IHshows the model predictions for these new values 
of f c . With this new effective cooling factor both WDM 
models (0.5 and 0.75 keV) provide equally good fit to the 
data at z = as the CDM model. On the other hand 
they now show the same over-prediction of galaxies at 



z = 1.5 as CDM. This is due to the very similar evolution 
of WDM and CDM models from z = 1.5 to z = 1 (sec 
Fig. 1 in Weinmann et al. 2012). 

Overall we find that by simply tweaking one single pa- 
rameter we have been able to get the exact same predic- 
tions in WDM and CDM models. This shows the danger 
in using a single observable (e.g. the stellar mass func- 
tion) to claim the success of one dark matter model over 
another. In the next section we will try to combine dif- 
ferent obscrvables to break the degeneracy between the 
DM properties and the parameterization of baryonic ef- 
fects in SAMs. 

3.2. Constraining WDM with the stellar-to-halo mass 

ratio 

In Fig. [2]we have shown that better match to the local 
stellar mass function can be achieved in the WDM mod- 
els by increasing the f c parameter in low-mass haloes. Of 
course this has the consequence of increasing the stellar 
mass content in low mass haloes. 

This increased ratio between stellar and halo mass can 
be constrained using several independent observations. 
The halo mass can be directly measured for stacks of 
galaxies using galaxy-galaxy weak lensing (e.g., Mandel- 
baum et al. 2006; Lcauthaud et al. 2012a, b) or satellite 
kinematics (e.g., Conroy et al. 2007; More et al. 2011), 
and inferred assuming halo abundance matching (Mostcr 
et al. 2010). A closely related quantity to the halo mass 
is the maximum circular velocity of the dark matter halo. 
The best observational probe of this is the rotation ve- 
locity in the "flat" part of the rotation curve, Vhat, as 
traced by 21cm observations of neutral hydrogen. We 
refer to the relation between stellar mass and outer ro- 
tation velocity as the Tully-Fisher relation, even though 
the original relation from Tully & Fisher (1977) was be- 
tween B-band luminosity and 21cm linewidth. 

In Fig. [3] we show the halo mass to stellar mass rela- 
tion at z = for our WDM models (with different / c ) 
and the CDM model. Data points with error-bars are 



■5 




i i i i i . i- 1 1 1 1 1 i i i 

10 8 10 9 10 10 10 11 10 12 

M star [M ] 

Fig. 3. — The stellar mass to halo mass relation. The data 
points are observational results using different approaches to de- 
termine the halo mass: Galaxy-galaxy weak lensing from Mandel- 
baum et al. 2006 (filled circles); Satellite kinematics from Conroy 
ct al. 2007 (open squares) and More et al. 2010 (open triangles). 
The lines are the predictions from the SAMs (CDM - solid black; 
WDM 0.75 KeV - red long-dashed; WDM 0.5 KeV - blue short- 
dashed). The dotted line shows the maximum stellar mass in given 
halo mass assuming an universal baryonic fraction. 

measured stellar to halo mass ratio from weak lensing 
and satellite kinematics (see Leauthaud et al. 2012a for 
a discussion of the different data sets.). The dotted line 
shows the maximum stellar mass in given halo mass as- 
suming the universal baryon fraction, /b = f2b /^m from 
the WMAP7 cosmology. 

While the three models give the same stellar mass 
function (as shown in Fig. [2]), they do predict differ- 
ent halo-to-stellar mass ratios especially for halo masses 
below 10 12 Mq. For a give halo mass the WDM models 
predict a large stellar mass with this difference increas- 
ing for decreasing m„ values. Unfortunately, the current 
weak lensing and satellite kinematics observations do not 
probe halo masses below 10 12 M Q . 

3.3. Constraining WDM with the Tully-Fisher relation 

For haloes with masses in the range 10 11 < 
^hoo/MQ < 10 12 , the best current probe of the halo 
masses comes from the Tully-Fisher relation. This con- 
straint is shown in Fig. |4j As before the observations 
are shown with points and error bars, while the models 
are shown with lines. For the observations the points 
show mean of log ilf sta r in bins of outer rotation veloc- 
ity, with the error bar corresponding to the error on the 
mean. There are typically ~ 10 galaxies per velocity 
bin. The observations are from the compilations of Stark 
et al. (2009) and McGaugh (2012). We calculate stellar 
masses using relations from Bell et al. (2003) (with a 
—0.1 dex correction to convert to a Chabrier 2003 IMF). 
We also convert the Hubble parameter from 75 to 70. 

For the models we compute the circular velocity at a 
radius enclosing 80% of the cold gas, i?so- At this radius 
model rotation curves tend to be flat (Dutton 2009). The 
SAM does not provide radial information for the stars, 



gas, or dark matter, so we use empirical constraints (from 
Dutton et al. 2011) for the stars and theoretical con- 
straints (Maccio et al. 2008; Schneider et al. 2012) for 
the dark matter. We assume the stars and cold gas are 
in exponential disks, with the scale length of the cold 
gas being 1.55 times larger than that of the stellar disk 
(Dutton et al. 2011). We determine the average disk size 
using the size vs stellar mass relation for late- type galax- 
ies from Dutton et al. (2011). Thus for a given stellar 
mass in the model, we know the average radius we should 
be measuring the rotation velocity, and the contribution 
of the stars and cold gas to this velocity. For the dark 
matter and hot gas we assume the profiles are NFW. For 
the CDM case we adopt the concentration mass rela- 
tion from Maccio et al. (2008) for a WMAP5 cosmology 
(which is very similar to that of the cosmology adopted 
here). For the WDM models we scale the concentration 
parameters according to the fitting formula of Schneider 
et al. (2012). At the radii we measure rotation velocities 
the galaxies are dark matter dominated, and thus halo 
response (whether it be contraction or expansion) has a 
minor impact on the predicted rotation velocity. 

Fig.|4]shows that the CDM model reproduces the data 
well, whereas at a fixed velocity the WDM models have 
higher stellar masses, especially at low velocities. At a 
velocity of Va a t ~ 80kms _1 the differences arc substan- 
tial: a factor of « 2 between CDM and 0.75 keV WDM, 
and a factor of « 1.6 between 0.5 and 0.75 keV WDM. 
These large differences are due to two effects. Firstly, as 
shown in Fig. [3] at fixed halo mass WDM models have 
higher stellar masses. Secondly, WDM models have lower 
concentration haloes, which results in lower rotation ve- 
locities, and hence higher stellar mass at fixed velocity. 

In summary the CDM model provides a good match 
to the data (note that the model was not tuned to 
match the TF relation), whereas both both WDM mod- 
els overprcdict the stellar masses at fixed velocity. The 
limit of m„ > 0.75 keV that we find here is consistent 
with current limits from large scale structure (e.g. Viol 
et al. 2008). 

4. CONCLUSIONS 

Recent observational results have challenged the other- 
wise successful Cold Dark Matter model on small scales. 
For example the inner density profile of dark matter halo 
is too concentrated to match the kinematics of satellites 
around the Milky Way (e.g. Boylan-Kolchin et al. 2011). 
Warm Dark Matter models have been suggested as a pos- 
sible solution since the suppression of power on small 
scales can in principle alleviate if not solve these issues 
(e.g. Lovell et al. 2012). 

However, all the CDM predictions are based on grav- 
itational only (N-body) simulations, which by construc- 
tion neglect the effects of baryons and their complicated 
network and interactions through gas cooling, star forma- 
tion and feedback that could possibly alter the results of 
pure DM simulations (e.g. Governato et al. 2012, Brooks 
et al. 2012) It is then important to investigate the pre- 
diction of CDM and WDM on the statistical properties 
of galaxy population, moving beyond a simple DM sim- 
ulation. 

In this paper, we performed high-resolution N-body 
cosmological simulations for three WDM models with 
m v = 0.5, 0.7 and 2.0 keV respectively, plus a controlled 
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Fig. 4. — The Tully-Fisher relation. The lines show the same 
models as in Figs. [Jfjj] E° r the models we calculate Vfl a t at a 
radius enclosing 80% of the cold gas. The filled circles with error 
bars show the mean (and error on the mean) of logM s tar in bins 
of Vfl a t using observations from Stark et al. (2009) and McGaugh 
(2012). The CDM model is clearly favored over both WDM models. 

CDM model. We couple these simulations with a Semi 
Analytical Model of galaxy formation to study the im- 
pact of a cut-off in the power spectrum on observable 
quantities such as the stellar mass function. 

For a fixed set of parameters describing the baryonic 
physics, models with low masses for the warm particle 
(m v = 0.5 and 0.75 keV), predict less galaxies with stellar 
mass Mgtar < 10 10 M Q than the current data at z = 
0. We find that the WDM model with m v = 2 keV 
provides almost identical results as CDM and it is able to 
successfully reproduce the data. The situation is reversed 
at higher redshift z = 1.5, with WDM models in better 
agreement with the observed stellar mass function than 
the CDM model (or the WDM model with m„ = 2 keV). 

However we show these differences in the stellar mass 
function, are strongly degenerate with the set of param- 
eters used in the Semi Analytical Models to describe the 
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(largely unknown) galaxy formation processes. By ad- 
justing a single parameter (the cooling rate of low mass 
haloes) we have been able to get indistinguishable results 
from all the DM models, both at high and low redshift. 
This shows that a single observable (e.g. the stellar mass 
function) can not constrain the effects of warm compo- 
nent on the galaxy formation process. 

In order to break this degeneracy we use independent 
constraints on the (integrated) star formation efficiency 
at low masses. We found that if the WDM models are 
tuned to reproduce the present (z = 0) stellar mass func- 
tion, the stellar mass for a given halo mass is system- 
atically larger in haloes of masses below (~ 10 12 M Q ). 
Such scales are not reliably probed by direct methods 
such as satellite kinematics or weak galaxy-galaxy tens- 
ing. In order to probe such low halo masses we use the 
relation between stellar mass and rotation velocity at 
large galactic radii, more commonly known as the Tully- 
Fisher relation. Current data already rule out models 
with m„ < 0.75 keV which is in agreement with other 
limits from large scale structure (e.g. Viel et al. 2008). 

Finally our study shows that by combining measure- 
ments of galaxy stellar mass function and stellar mass 
- halo mass relation down to low galaxy masses (« 
70km/s) it is possible to obtain very tight constraints 
on the mass of a possible warm component. This opens 
a new window in the search for the nature of the elusive 
dark matter component of our Universe. 
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